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Abstract: 

There exists experimental evidence that a dibaryon resonance d' with quantum num¬ 
bers J^=0“, T=0 and mass 2065 MeV could be the origin of the narrow peak in 
the {n^,7r~) double charge exchange cross-sections on nuclei. We investigate the 
six-quark system with these quantum-numbers within the constituent quark model, 
with linear conhnement, effective one-gluon exchange at short range and chiral in¬ 
teractions between quarks (vr- and a-exchange). We classify all possible six quark 
states with J^=0“, T=0, and with N=1 and N=3 harmonic oscillator excitations, 
using different reduction chains. The six-quark Hamiltonian is diagonalized in the 
basis including the unique N=1 state and the 10 most important states from the 
N=3 shell. We hnd, that with most of the possible sets of parameters, the mass of 
such a ’’dibaryon” lies above the N(939)-|-N*(1535) threshold. The only possibility 
to describe the supposed d’(2065) in the present context is to reduce the conhne¬ 
ment strength to very small values, however at the expense of describing the negative 
parity resonances N*. We also analyze the J^=0“, T=2, N=1 six-quark state. 


* Alexander von Humboldt Fellow, on leave from Alma-Ata Power Engineering In¬ 
stitute, Kosmonavtov 126, Alma-Ata, Kazakhstan. 


1 Introduction 


Believing QCD to be the theory of strong interactions, we consider hadronic systems 
as built up by quarks and gluons. With the exception of lattice calculations, which 
start directly from the underlying current quark and gluon helds, the descriptions 
of baryons make use of some effective degrees of freedom of nonperturbative QCD, 
such as constituent quarks |^, chiral helds 0, etc. In many-baryon systems, like 
atomic nuclei, the relevant degrees of freedom are baryons and mesons (except for 
perhaps some special cases 0, where the quark degrees of freedom must be taken 
into account explicitly). 

However, in principle nothing forbids the existence of objects with baryon number 
bigger than one, which have some QCD motivated origin and could be considered 
as built up with quarks and gluons rather than baryons and mesons. Such objects 
would have, in contrast to atomic nuclei, quite a small size, of order of 1 fm. 

This idea was very popular about 10-15 years ago and a lot of efforts were 
made to hnd such objects, both theoretically and experimentally (see, for example, 
]I0[| and references therein). However, an experimental search for dibaryons 


refs. I 


in the NN-system has not been successful until now. 

The reason for this seems to be quite clear. The only demand for the color part 
of the six-quark wave function of a hypothetical dibaryon is that it is a color singlet, 
i.e. it has a symmetry 

|2"Ig = ffl ■ (1) 

I I IG 

This singlet SU{3)c representation contains in its Clebsch-Gordan expansion the 
state 

[l']c ^ [l"]c = B 0 B (2) 

□ c Dc 

(Here and in the following, ”x” denotes an inner product, whereas ”0” denotes 
an outer one). But in this case the strong conhnement forces ~ A“ ■ A“ are absent 
between two color singlet objects: With the usual normalization of the Gell-Mann 
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matrices = A“A“ = ^ (summation over double indices), one finds 
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) = 0 
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*C{1,2,3} 
jC {4,5,6} 


(3) 


So, if two baryons are very close to each other and form a 6q-system, the con¬ 
finement forces will not prevent such a system from a decay into baryon-baryon 
channels. This means, that the lifetime for such a system would be extremely small 
and the corresponding large width would thus not allow to treat such a system as 
a dibaryon. The only possibility for the dibaryon in this case would be if there 
were some forces (like A“ • A“<Tj ■ (Tj , r* • TjCTi ■ (Tj, etc.) which bind the six-quark 
system below the corresponding baryon-baryon threshold for a strong decay. This 
argument stimulated for example the search for an H-particle i a M- 

On the other hand, there seems to be some recent experimental evidence for a 
negative parity resonance in the vrA^A^-system with the quantum numbers J'^=0“, 


T=0 1^, for which the notation d’ has been introduced. The pionic double charge 
exchange reaction (DCX) on nuclei 


(7r+,7r-) ^Z + 2 


( 4 ) 


exhibits, independently of the nuclear target, a very narrow peak near the pion 
kinetic energy T^r = 50 MeV. Its position as well as its width are practically identical 
for the available world DCX-data on light and medium nuclei such as 
'^"^Ca, ^®Ca. Only its amplitude depends on the considered nucleus. This suggests, 
that it can only be connected to some elementary process. 

Due to charge conservation this elementary process involves at least two nucleons 
within the nucleus. The hypothesis, that this peak is due to a dibaryon resonance d' 
in the nNN channel with quantum numbers J^=0“, T=0 and mass 2065 MeV [pT 


allows to describe all available data. This dibaryon cannot decay into the nucleon- 
nucleon channel due to the Pauli principle and its mass is below any baryon-baryon 
[NN*, N*N *,...) threshold. These peculiarities explain the very small width of only 
T ~ 5 MeV of the d’ dibaryon. 
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These circumstances motivated our detailed study of the six-quark system with 
the quantum numbers J^=0“, T=0 within a constituent quark model P, |^, |T^, 


151 with and without chiral interactions between constituent quarks. 

Dibaryons with quantum numbers J^=0“, T=0 were already studied in the 
framework of string-like (deformed) bag models and were calculated as strings with 
and colored quark clusters at the ends [|, An essential shortcoming of 
this model is the lack of antisymmetrization between quarks belonging to differ¬ 
ent clusters. Antisymmetrization can be neglected only for well separated clusters. 
However, in the — system the color-exchange (conhning) forces between quarks 
do not allow large separations, and the typical size of such a system is expected 
to be of order 1 fm. We know from experience in the A^A^-system that the Pauli 
principle on the quark level plays a decisive role at such distances f 
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The other drawback of the bag-model in the six-quark system is that the bound¬ 
ary conditions (which simulate the quark confinement) prevent such a system from 
the color-singlet 3q - color-singlet 3q clusterization, which is the most important 
phenomenon in 3n-quark (n > 1) systems. If such a model were correct, a large 
number of dibaryons should be observed. From this point of view, a potential model, 
where quark conhnement is approximated by two-body q-q forces ~ A“ ■ A“, is more 
satisfactory, although one has to take care of the long-range Van-der-Waals forces 
in this case, when describing the NN-system . 


In this work we study the six-quark system with the quantum numbers J'^=0 , 
T=0 within the constituent quark model i 0, [I^, |T^ taking into account 

properly the Pauli principle throughout the whole calculation. We investigate both 
possibilities; with a chiral held (pions and sigma mesons) coupled directly to con¬ 
stituent quarks, and without such a chiral held. In the latter case, quarks interact 
only through color-exchange potentials. In our previous communication |^, only 
the lowest shell-model state s^p was used. In that paper, we did not care about a cor¬ 
rect description of the one-quantum nucleon excitations N*(1535) (J^=^ , T=|). 
However, this is quite important, since the dibaryon mass should be compared with 
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the N(939)+N*(1535) threshold. Here we present a more complete calculation tak¬ 
ing into account conhguration mixing, including the 10 presumably most important 
excited states. We hnd that with all sets of parameters, that describe rather well the 
lowest baryon mass spectrum, the mass of the J'^=0“, T=0 6q-state lies around 100 
- 200 MeV above the NN*-threshold. However, including more excited states in our 
basis, this mass could come down below the threshold by some 10 MeV’s. If it is so, 
there really could exist a dibaryon with these quantum numbers. On the other hand, 
its mass, at least within our model, is still higher than the experimentally observed 
peak suggests. Assuming that the conhnement potential in a six quark system is 
weaker than the usual choice in three quark systems (This would still describe the 
nucleon and the A, but not the negative parity N* resonances.), yields a dibaryon 
with a mass close to 2065 MeV and an oscillator length fcg — 1-25 fm. 

It seems that this 6q-system is very sensitive to the explicit form of the con¬ 
hnement mechanism, or, in other words, the extension of the baryon conhnement to 
multi-quark conhnement might not be too straightforward. In this sense, the con- 
hrmation of the d’ resonance could give valuable information to our understanding 
of conhnement. 

The structure of the paper is the following: The ehective Hamiltonian for the 
dibaryon is presented in section 2, and some short motivation will be given for its 
diherent ingredients. Fitting the parameters to the spectrum of light baryons is 
the commonly accepted procedure to hx the ehective Hamiltonian. In section 3, 
we present the Translationally Invariant Shell-Model (TISM) basis. We classify the 
basis states for the 6q-system with J^=0“, T=0, and review some important points 
concerning the fractional parentage technique for the wavefunction. We present 


the results for the dibaryon, hrst excluding the mixing of excited states (cf. |[T7||), 
and then including the ten presumably most important excited states. A short 
discussion of a possible J^=0“, T=2 six-quark state is also added. Section 7 gives 
a conclusion of the present work. All necessary analytical expressions are reserved 
for the appendices. 
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2 The effective Hamiltonian 


One now generally assumes, that the constituent quarks are quasiparticles with a 
complicated structure and dynamical (q-dependent) mass The current 

quarks of the underlying QCD acquire the constituent mass due to the spontaneous 
breaking of the chiral symmetry of the QCD Lagrangian, which is nearly exact in 
the SU (2)-flavour sector. 

This chiral symmetry breaking is caused probably by the instanton structure 

, |2^ and is characterized by the corresponding order pa- 


of the QCD vacuum 


rameter - the non-zero quark condensate - and by the Goldstone exitations - the 
pions. All these features are also reproduced in the Nambu-Jona-Lasinio model (see 
recent review and references therein), which can be considered as an effective 
approximation to the underlying QCD in the low-energy domain. 

So, the low-momentum Lagrangian must contain the constituent quarks and 
chiral helds as effective degrees of freedom. But it must also contain a color-exchange 
interaction between the constituent quarks, since the latter are coloured objects. 
Among these effective interactions we should include a conhnement force (whose 
real nature is not known), and also a short-range (i.e. at distances less than the 
chiral symmetry breaking scale) effective gluon-exchange between quarks. 

There is a very close analogy to this picture coming from solid state physics. The 
correct effective degrees of freedom for an explanation of thermo- and/or electric 
properties of metals are not light elementary electrons and ions, but heavy electrons 
with effective mass m* and phonons. Phonons are Goldstone excitations and ap¬ 
pear as a result of the breaking of translational invariance in the lattice of ions. This 
spontaneous breaking is also characterized by the effective electron mass m* (micro¬ 
scopically, this effective mass arises due to very complicated interactions between 
the elementary electrons and the lattice). There are also some residual interactions 
such as electron-electron, electron-phonon, etc. 

After these quite general remarks and motivation for our chosen quark picture. 
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let us now turn to the concrete description of the effective Hamiltonian. 


The chiral sector: Pion— and Sigma—exchange 

The chiral-invariant interaction Lagrangian of the cr-model type [^, ^ p!4| , 
which describes in our case an interaction of the constituent quark field T (spinor 
in SU{2)fiavour^space) and the chiral helds tt and a, is 

Ant = -fi'd'(cr + zysTT • t)T. (5) 


The constituent quarks have a complicated structure which can be parametrized 
by a formfactor in the quark-chiral held vertex. For that, we use the substitution 


9 


9 


( 6 ) 


where A characterizes the chiral symmetry spontaneous breaking scale, and must be 


of the order of the instanton size in the QCD vacuum |^, i.e. A ~ (0.7 — 1.5)GeV 
This cut-off contributes to the pion-nucleon interaction radius as 

3 


{n.N = b^ + 


A2 


(7) 


where b is the quark core size of the nucleon [|T^. The larger the cut-off A, the 
more pointlike behaves the pion-quark coupling. Intuitively, one could think of this 
cut-off as an effective size of the constituent quarks (r^) = For example A=4.2 
fm“^ would correspond to a size of \J— 0.4 fm. On the other hand, the cut-off 
formfactor behaviour in eq.(^ can be related to the g^-dependence of the dynamical 


constituent quark mass [14 


The sigma mass ruo- is fixed by the relation rria ~ 2mg, where rUg ~ \mN is the 
constituent quark mass. On the other side, the relation 


2 /t 2 , 2 

= Arrig + 


appears as a result of the bosonization of the Nambu and Jona-Lasinio Lagrangian. 

The pion-quark coupling constant g is determined by the well known cou¬ 
pling constant g = 9 -kN {9nN = 13.36), where the coefficient | comes from the 
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spin-isospin matrix element when we consider the vriV interaction as the interaction 
between the pion and the 3 constituent quarks. 

Our notations for the coupling constants are 


~ 2m 


9crq 


mo 

2m, 


■9-9 


uq A.,,uq 

Fourier-transforming the static pion propagator, 

1 f , ^ ^ 1 exp(—pr) 

- I exp{iq ■ r) ^ ^ =- 


( 8 ) 


( 9 ) 


(27r)3 7 “ ^ ^2 _j_ ^ 4^ ^ 

we arrive at the usual Yukawa-like potential for the isovector pion (and similarity 
for the scalar-isoscalar sigma, which is not sensitive to any spin-isospin quantum- 
numbers of the quarks, but only to their orbital distribution) : 


= 




Ti ■ Ti 


nq ' i 


h? — m? dvr 

TT 

A2 


o-i ■ o-jVc{m^r) + SijVT{m^r)) 


A2 


mt 




► A) 


= ~ A2_^2 ^{Vcirriar) - Vb(Ar)) , 


Vc{mr) = 


exp(—mr) 


r 


Vrirnr) = [l -h — H- 

V mr m^r^ / 


exp(—mr) 


_ , 30 -,-r (Tyr 

I ^2 


( 10 ) 

( 11 ) 

( 12 ) 

(13) 


The Color—part: Confinement and Gluon exchange 

We take the conhnement potential, which mainly determines the medium- and 
long-range (on the quark scale) phenomena, to be linear 

(14) 


r.f'"-'(r=|r,-ri|) = -Yty(a,r-C) 


The effective ” one-gluon exchange” potential, which is responsible for the very 
short-range phenomena and motivated in its form from the charmonium spec¬ 
troscopy 1^, contains essentially hve terms: 




OGE 




■S',, 




r X pij ■ (cr, cTj) 


(15) 
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the so-called color-coulomb (CC ~ ■ A“ A), the contact color-electric (CE ~ 

Af ■ Aj S(fJ) term, the contact color-magnetic (CM ~ A“ ■ A“crj ■ ctj S{f)) term, the 
tensor term proportional Sij and the Galilei-invariant spin-orbit term. 

In our present study, we drop the spin-orbit term. As it is well known, there 
is practically no room for the spin-orbit interaction when describing the baryon 
spectrum. It is well seen, for example, from the small mass-splitting of the A 
and I resonances N*(1535) and N*(1520). Usually it is assumed, that the Galilei- 
invariant spin-orbit forces from the OGE-potential are compensated in the baryons 
by the spin-orbit forces from other sources, e.g. the Thomas term coming from the 
conhnement [Q or the cr-exchange 

Parameter fitting to light baryons 

The constituent quarks, tt- and cr-helds, and the color-exchange forces are some 
effective degrees of freedom of nonperturbative QGD. They do not follow from hrst 
principles in a compelling and unique way, but are only supposed to simulate the 
dominant features of QGD in the low energy domain. The strong coupling constant 
as for example, is not the momentum-dependent running coupling constant of QGD, 
but an effective parameter, which has to be hxed to some observable of low energy 
QGD to define our Hamiltonian. 

A possible way to £t the 5 parameters of our final Hamiltonian of eq. ( 0 ). 
as,ac,C,A and mg, is to describe the baryon spectrum, or more precisely the light 
baryons such as the nucleon, the A-resonance, and the J^=| )T=A one-quantum 
excitation N*(1535) (the reason, why one should fit our parameters first of all to 
the N*(1535), but not to the other excited states, will be seen from considerations 
in section 4). For the N and A, we use the simple harmonic oscillator functions 
(without the center-of-mass oscillations: see states \Nq) and |Ao) in appendix A 
and in refs.^, [^). The analytical expressions for the mass-formulae mjsi = and 
can be found in appendix A. 

To suppress the mixing with the two- and more-quantum harmonic oscillator 
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states in the N and the A, we impose the nucleon stability condition []^ = 

0, where is the harmonic oscillator parameter in the 3-quark system (which 
coincides with the mean-square quark core radius of the N and the A). 

For the N*(1535), we diagonalize analytically our Hamiltonian in the space in¬ 
cluding the two possible negative parity one-quantum states and in 

notations of the second paper in ref.p| (All the necessary wave functions can be 
found there), and £t the smaller eigenvalue to the N*(1535) mass. The correspond¬ 
ing mass-matrix can also be found in appendix 7\Q. 

Thus, for the six free parameters, as,ac,C,A,hN,mq, we have now four con¬ 
straints. But there is not to much freedom because one has the following additional, 
approximate constraints: 

i) the parameter A is hxed by the scale of the chiral symmetry spontaneous breaking, 
A ~ 0.7 - 1.5 GeV, 

ii) the quark core radius of the nucleon should be in the region — 0.45 - 0.6 fm, 

iii) the constituent quark mass should be rriq ~ 200 - 400 MeV. 

In table |^, we show possible sets of parameters. The sets I-III correspond to 
the inclusion of all types of the q-q forces, while in set IV, there are no chiral helds 
inside the baryons, i.e. no tt- and a-exchange between quarks. This selection does 
not give a systematic study of how the results depend on the chosen parameters, 
but is meant to show some qualitative features of our results. 

Increasing up to bjsi = 0.6 fm in the above discussed parameter htting scheme, 
we get a too small constituent quark mass, ruq ~ 170MeV. 

We also show in this table the model mass for the two-quantum positive parity 
J^=|^,T=i resonance. This mass is calculated by diagonalisation of our Hamilto¬ 
nian with already hxed parameters in the space including two possible two-quantum 
excitations |Ai) and IA 2 ) with L=0 (see ref. 0 and appendix A). 

We see that this mass is essentially higher than the mass of the Roper resonance. 


^In this basis, the 2nd eigenvalue can be compared to the next one-quantum excitation with 
J^=i ,T=i, the N*(1650). We give this result for illustration in table ||. 
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The Roper resonance is a longstanding problem in the description of baryons within 
the constituent quark model. We would like to remark at this point, that there are 
some new ideas for the CQM, that seem to overcome this problem. 

Set V finally gives a set of parameters, determined from the following fitting 
scheme, which will give the dibaryon results closest to the ’’experimental value”. To 
obtain these parameters, we fix as input the constituent quark mass niq, the cut-off 
A and the nucleon hadronic size to determine as, ac, C by the A-N mass splitting, 
the nucleon stability condition, and the nucleon mass respectively. In the table we 
see the small conhnement strength Oc — 25 MeV/fm, which goes together with a 
bigger quark core size of the nucleons 6Ar=0.6 fm. These parameters emerge also 
naturally, if the nucleon quark core size is determined by the nucleon stability 
condition. A similar parameter set is for example used to describe the nucleon 


magnetic moments in |T^. However, it is not possible to simultaneously describe 
the one-quantum excitations with these parameters. 


So, with fixed parameters, the effective 6-quark Hamiltonian is given by 


//'«’>=u, + 


z=l 


2m, 


p2 6 

* cm I \ ' -irlnt 

12m ^ ^ 

i<j 


(16) 


where 


ylnt _ yConf yOGE _|_ y-K _|_ yo 


(17) 


This Hamiltonian contains interactions of two particles only, so the two-particle 
fractional parentage coefficients (fpc), which will be reviewed below, are particularly 
well suited to evaluate matrix elements for two-body interactions of the above type. 


3 The Translationally Invariant Shell Model basis 
(TISM) 

The harmonic-oscillator basis is highly successful in describing the baryon spectrum. 
Let us recall, why this is so. The reason is, of course, that the quark-quark harmonic 
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oscillator attraction simulates well the quark-quark confinement forces in an 3-quark 
color-singlet object. Indeed, the color-singlet 3-quark wave function 


[lie = 


c 


is antisymmetric in color space in any quark pair. Since 

8 


(Sci^“ ■ Sd 


(18) 


(19) 


the confinement potential of eq.(|^ is attractive in any quark-quark pair. Of course, 
historically, the confinement potential was just constructed as to provide an attrac¬ 
tion. 

In the six-quark system, the situation is very different. Here the color part is 
characterized by the mixed permutational symmetry 


[2lc ^ 


c 


( 20 ) 


It means, that there are both antisymmetrical and symmetrical pairs. But in the 
symmetrical ones, instead of attraction, the potential (|^ leads to repulsion, since 

4 


{[III]c|A“-A“| [IHc)^ 


( 21 ) 


So, the confinement forces (H) try to build the six-quark system in such a way as 
to divide it into two color-singlet 3-quark clusters. 

It is not possible to simulate this behaviour with only one harmonic oscillator six- 
quark configuration. In the harmonic oscillator wave function we have an attraction 
in each quark-quark pair|. However, one can use the 6q harmonic oscillator as a basis 
to diagonalize the Hamiltonian eqs. (ED,®)- In this case, it is intuitively clear, that 
the confinement force will mix the different 6q-harmonic oscillator states in such a 
way as to provide in the final six-quark wave function a color-singlet - color-singlet 
clusterization. If the baryon-baryon asymptotics with well separated baryons (e.g. 
for the deuteron) exists in some six-quark system, the harmonic oscillator basis is 
^Due to the same reason, the bag model for a six-quark system is conceptually not correct. 
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not efficient to describe this system at all range, since we would need a tremendous 
number of highly excited harmonic oscillator configurations in this case[|. 

In the 6q-system with quantum numbers J^=0“, T=0, the nucleon-nucleon 
clusterization is suppressed due to the Pauli principle on the nucleon level. On the 
other hand, if there is a deeply bound state with respect to the N(939)+N*(1535) 
threshold, the harmonic oscillator basis could be quite successful in describing such 
a system. 

The shell-model basis allows, as we will see in more detail in the following, to 
respect the Pauli principle, i.e. the antisymmetrization of the total wavefunction, 
at any level of the calculation, so that effects arising due to the Pauli-principle are 
naturally included in the formalism. This is very important for interquark distances 
around or smaller than 1 fm, which are expected for a dibaryon. The importance 
of the Pauli principle is well known from the NN-system. There, the short range 
repulsion at distances smaller 1 fm is essentially due to quark exchanges, caused by 


the Pauli principle on the quark level |12 


Furthermore, another advantage of this Hamiltonian lies in the fact, that the 
center-of-mass (cm) motion is removed properly, so as to exclude spurious states 
from the harmonic oscillator basis. 

The price we pay for the exact inclusion of the Pauli principle and unambigious 
exclusion of spurious states in our basis, is of course the uncertainties induced by the 
non-relativistic description and by the effective nature of our Hamiltonian, which is 
not derived from first principles, but is assumed to simulate the dominant features 
of low-energy QCD. 


We remind the reader that the harmonic oscillator six-particle states are exact 

^But if the six-quark harmonic oscillator states at small distances are combined with the cor¬ 
responding clusterized configuration like 2, 3)52(4, 5,6) x(r)} at medium and large range 

in one variational task, we do not need much harmonic oscillator states at small distances. Such a 
basis is essentially more flexible than the one-channel RGM basis. This variational program was 
realized in part in the investigation of the NN problem in refs. MM- 
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solutions of the Translationally Invariant Shell Model (TISM) Hamiltonian, ref. 



( 22 ) 


The harmonic oscillator-eigenstates of this Hamiltonian 


|h.o. — state) = \A = 6 , N [f]x (A/r), L, S,T, a) 


(23) 


are classihed (cf. ref. [^) by the following set of quantum numbers; N indicates 
the number of internal excitation quanta, the Young-pattern [f]x determines the 
spatial permutational symmetry, the Elliot symbol (A/r) gives the SU{3) harmonic 
oscillator multiplet, L, S,T are the total orbital angular momentum, total spin and 
total isospin. The total spin S is uniquely connected with the Young pattern [f]s 
for the SU{2)s representation by S' = (/i — / 2)/2 where /i and /2 are the hrst and 
second rows in the Young pattern [f]s- The same applies to the total isospin T. 
The unambigious dehnition of a h.o.-state requires, as will be shown in a moment, 
additional quantum numbers, here collectively denoted by a. 

In conventional nuclear physics (dealing with nucleons), due to the Pauli princi¬ 
ple, the spatial Young pattern [f]x must be conjugate to the SU{4:)st Young pattern 
[f]sT in order to provide the total antisymmetry (i.e. [f]x = [f]sT, where the tilde 
implies the Young pattern reversed with respect to the main diagonal). 

In the quark model (with two flavours), due to the additional degree of free¬ 
dom color, the classihcation is more complicated. For a six-quark system, the 
full permutational color symmetry is described by the [2^]c color Young pattern 
{SU{3)c color singlet). To identify unambigiously the six-quark state we have to 
use additional quantum numbers, because the inner product of the SU{3)c color- 
singlet representation and the SU{2)s spin representation (or the SU{2)t) con¬ 
tains in general more than one representation of the SU{6)cs group 


SU{6)cs 2) SU{3)c X SU{2)s (or SU{6)ct)- The spatial symmetry [f]x determines 
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uniquely only the SU{12)cst permutational symmetry, [f]x = [/]csr; so we need 
in addition [f]cs or [f]cT- Other reduction chains are also possible, for example 


with the intermediate spin-isospin SU{A)st symmetry |^, or with the intermedi¬ 
ate spin-spatial symmetry [^, etc. All these reduction chains are equivalent, and 
the choice is determined by convenience (for different types of quark-quark forces 
different reduction chains are convenient). 

In our preliminary communication |]^, we have used the SU{6)cs and SU{4:)st 
symmetries to classify all possible states. In the present paper the choice of the 
intermediate spin-spatial permutational symmetry [f]xs is entirely motivated by 


the available tables of two-particle fractional parentage coefficients in refs. [31 


Classification of the states with J^=0 , T=0. 

We get the negative parity only if our states contain an odd number of harmonic 
oscillator excitations, i.e. the possible number of excitation quanta N is restricted 
to 

TV = 1,3,5,... (24) 

There is only one state with N=l, compatible with J^=0“, T=0: 

|iV = 1, [51]x, (A/i) = (10), L = 1,S=1,T = 0, [42]xs) • (25) 

In this case, all possible intermediate permutational symmetries (within different 
reduction chains) are determined simultaneously and uniquely 

[321]c5 , [321]sr , [2211]ct , [321]xr , [42]xs , [321]xc • (26) 

This harmonic-oscillator state (i.e. TISM-state) is uniquely connected with the 
harmonic-oscillator shell model state s^p as follows: 

\s^p [51]x (10)L = 1) = $ooo(Rcm) ■ |iv = 1 [51]x (10)L = 1) (27) 

The shell model state s^p with the spatial symmetry [6]x is a spurious one, since 

\s^p [6]x (10)L = 1) = <l>iii(R,^) ■ |iV = 0 [6]x (00)L = 0) (28) 
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There are a lot of different states with three internal excitation quanta. 

The full classihcation of the orbital parts with N=3 and the spatial symmetries 
[42]x, [411]x, [33]x, [321]x, and [3111]x was done in Ref. [p7[] . In addition we have 
also the states with [51]x spatial symmetries (both [6]x and [51]x are also allowed 
in the quark-model, but the state with N=3, [6]x is a spurious one). 

However, not all of these orbital states are compatible with the color, spin and 
isospin quantum numbers of the d'. In table we present all possible conhgurations 
with N=3, compatible with the quantum numbers J'^=0“, T=0. 

As mentioned above, all different reduction chains are equivalent. For example, 
for the two reduction chains 


SU{12)cst 3 SU{6)cs x SU{2)t D SU{3)c x SU{2)s x SU{2)t , (29) 
SU{12)cst 3 SU{3)c X SU{4)st 3 SU{3)c x SU{2)s x SU{2)t , (30) 


the intermediate SU{6)cs or SU{4:)st symmetries are needed respectively to define 
the shell-model states uniquely. They are given in line 4 and 5 of table ^ (see for 


more details ref. 17 


The reduction chain we hnally adopt in our calculation, decouples color and 
isospin from the remaining spin-orbital symmetry. 


[l^]xc5T [33 ]t X [f']xcs [33]t x [222]c x [f]xs 
D [33]t X [222]c X [f]s x [f]x 


(31) 


In the sixth line of table |^, we give the intermediate [f]xs symmetry, needed in this 
classihcation scheme. 

The rules for inner products of Young patterns for certain permutational sym¬ 
metries can lead to the situation, that a certain product symmetry occurs more 
than once. In these cases we need an additional quantum number (in table |^, these 
quantum numbers are indicated by subscripts) to distinguish these different states. 

The last line of table ^ gives the number of states in each column, all in all, we 
hnd 31 TISM-states with N=3 and the quantum numbers of the d’. Even if it were 
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possible to include all these 31 excited states in our calculation (this would give a 
32 X 32 matrix problem), this would not be a complete basis, since we are not taking 
higher excitations with N = 5,7,... into account. For a qualitative estimate of the 
effect of configuration mixing with the N=3 states, we choose the 10 states in the 
third row of table |^, characterized by 

N = 3,[42]x,L = l,S = l,T = 0 , (32) 


in order to build a basis of 11 states for our calculation. The notation for any of 
these ten states is analogous to eq. (^) and can be easily written down with the help 
of table |. 

The reason why we choose these states is, that there are the two most ’’symmetri¬ 
cal” (i.e. the longest) Young patterns [f]cs and [f]sT among these states. We remind 
the reader, that the more symmetrical pairs in color-spin-space we have, the more 
attraction in the six-quark system we get, due to the color-magnetic A“ ■ A“crj • a-j 
forces in eq. ( 0 ). 

This is well seen from the matrix elements 


{[fij]cx[fij]s\ \[fij]cx[fij]s) - { rni*’ (33) 

The same situation takes place for the pion-exchange forces (^). Here, one has 


4 

3’ 

if [2]c, 

[2]5, 

[2]cs 

8, 

if [ll]c, 


[2]cs 

-4, 

if [2]c, 


[ll]cs 

-8 

3 ’ 

if [ll]c, 

[2]5, 

[lljcs 


X [/p]t| T* • TjCr^ ■ (Tj I [fij]s X [fij]T) 


1, if [2]s, [2]r, [2]sr 

9, if [ll]s, [11 ]t, [2]sr . 04 ', 

-3, if[ll]5, [2 ]t, [11]5t 

—3, if [2]5, [11 ]t, [lljsT 


So, at quark-quark distances larger than, for example, 1 fm, where the first term 
in eq. (pT]) dominates over the third one, one has an attraction due to central tt- 
exchange forces in antisymmetrical FT-pairs. At small quark-quark distances (less 
than 1 fm), the cut-off in eq. (|IOD dominates, and one has an attraction in symmet¬ 
rical S'T-pairs. Since the characteristic size of our system is less than 1 fm, one has 
an attraction in symmetrical S'T-pairs. 
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The fpc (fractional parentage coefficients) technique 


Here, we cite briefly the definition and the application of the fractional parentage 
expansion (fpc), refering the reader to the corresponding literature |^-|^ for more 
details. 

Each six-quark harmonic oscillator configuration can be presented by means of 
the fractional parentage technique and the Talmi-Moshinsky transformation as a 
superposition of various four-quark x two-quark x relative motion components: 


\NaLSJT)= Y. ■.LSJT\ (35) 

Here, a stands for the necessary set of additional quantum numbers, including 
[f]x) (A/i), [f]xs, etc. to define unambigiously the six-quark state, and ^Bi,B 2 ,n,i 
is the so-called fractional parentage coefficient (fpc) in the TISM. 

(pBi and 0^2 are the TISM states for the first four particles and for the last two 
ones respectively. (pni{'^ is the n-quantum harmonic oscillator function with 

^ r 1 + f2 + r 3 + r 4 r 5 + r g 


The summation in eq. (p5D is carried out over all possible internal states of clusters 
Bi and B 2 and their relative motion n, I, provided that 


A^i + A^2 + ^ 


The specific feature of expansion (|^) is, that the antisymmetric six-quark wave 
function is expanded into the sum of orthogonal, but not fully antisymmetric terms. 
Each term is antisymmetric only within the clusters Bi and B 2 . 


By use of expansion (^), the interaction energy matrix elements of a pure two- 


body potential Hint = Vij can directly be evaluated: 


^i<j "o 


( N a LSJT\ Hint |iV a' L'S'JT) = x 


(36) 


X E 

Bl ,B2,B2 


r“ r" V 

Bi,B2,n,l ^ Bi,B2,n,l ^ 


( angular 
momentum 
recoupling 
V matrix 




X (0B2(5,6)|1/56|0b,(5,6)) 
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The expression (|36| ) is written in symbolical form to avoid bulky 6j- and 9j-symbols 
appearing due to the necessary recoupling of the angular momenta when calculating 
a contribution of the non-central forces. Here, a summation over all intermediate 
momenta is also assumed. For central forces, the recoupling matrix is absent. 


It will not be our task here to rederive the rules of how to construct the two- 
particle fpc’s for the 6-particle system. General considerations for the construction 
of fpc’s can be found in various publications, as for example in |^-|^ for the 
6-quark system, and the fpc’s (or more precisely the scalar factors) needed in our 
case are tabulated in ^ and |^. It is shown in these articles (and in the 

articles cited therein), that for a given reduction chain, the total fpc-coefficient 
factors out in a product of several scalar factors of the Clebsch-Gordan coefficients 
of the corresponding unitary groups, each one associated with a single step of the 
reduction. 

In our case, the corresponding fpc are given as a product of scalar factors 


_ Cr'U Q Q Q TP^ F 

^ v.Q ' or ^ ■ i 


XSC-T 


XS-C 


X-S 


C 


TISM 


(37) 


The hrst factor in eq. ( 0 ) is the weight factor and is determined only by the 
dimensions of the irreducible representations of the permutational group in isospin 
space for six particles (Young pattern [/]r) and for the hrst four particles in eq.(^^ 
(Young pattern [/i]r): 


SF 


u 

XSC-T 


IT 


\ dim[/]^ 


The next two factors in eq.(^) can be extracted from the tables 2 - 5 in ref. |pi 


The color scalar factor (i.e. the scalar part of the SU{3)c Glebsch-Gordan coefficient 
in the reduction SU{3)c Y 0{3)c) can be found in table 1 of ref.. 

Finally, the last factor in eq. (|37D Ttism, he. the orbital part of the fpc, can be 
found in table 2k of the hrst paper from |^. This table contains all the necessary 
coefficients for N=3, [f]x = [42]x six-particle harmonic oscillator conhgurations. 
For the lowest conhguration, i.e. for |Y = 1 [51]x (A/r) = (10) L = 1, S' = 1,T = 0), 
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these coefficients are trivially calculated, and one finds: 

{N = l[51]x{W)L = l\N, = 0[A]x{00)L, = 0;nl = 00,N2 = l[ll]xL 2 = l) = -1 
{N=l[51]x{l0)L = l\Ni = 0[4]x{00)L^ = 0-,nl = n,N 2 = 0[2]xL2 = 0) = -1 
(Ar = l[51]x(10)L = l|Ari = i[31]^(10)Li = l;n/= 00,Ar2 = 0[2]xh2 = 0) = 1 

(39) 


4 Results for the single N=1 configuration 


As our full six-quark Hamiltonian of eq-dlg) for the dibaryon comprises not only 
confinement, but also effective gluon-, pion- and sigma-exchange, it is very different 
from the harmonic oscillator Hamiltonian of eq.(E 


In our previous paper |]^, we treated the difference between the effective Hamil¬ 
tonian of eq. (p!6D and the harmonic-oscillator Hamiltonian of eq.(p^) as a residual 
interaction. The corresponding diagonal matrix element for the lowest possible con¬ 
figuration, given by eq.(p^) 

= = 1, [51]x(10)LSTJ = 1100|ff<“''>|iV = 1, |51]x(10)LSTJ = 1100) , 

(40) 

can be found in ref.[|^] and also in appendix B. 

As we already mentioned in the previous sections, the confinement forces of 


eq.(|T^ work towards the color-singlet -|- color-singlet clusterization in the 6q- 

system. In our case, that could be NN*_(1535) and NN3_(1520) clusterizations. 

2 2 

We can exclude higher excited nucleon states or N*N* clusterizations, since their 

threshold energies would be quite high. However, the NN*_(1535) is highly prefer- 

2 

able compared to the NN3_(1520). The reason is, that the latter cluster component 

2 

in the J^=0“, T=0 six-quark system is not compatible with the lowest possible 

harmonic oscillator configuration of eq. (P3|) . Indeed, in the NN3_(1520) system, 

2 

the relative motion angular momentum must be Z = 2 to provide total angular mo¬ 
mentum and parity J^=0“ of the six-quark state. Such an angular momentum I 
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needs a non-zero number of harmonic oscillator quanta, n = 2,4,..corresponding 
to the relative motion. So, together with the one quantum from the internal N* 
excitation, the total number of harmonic oscillator quanta would be at least 3. The 

NN* + (1440) component is also incompatible with the N=1 six-quark state. Thus, 

2 

the ’’dynamical” threshold for a possible J^=0“, T=0 dibaryon is N(939)+N*(1535). 

It is instructive to compare the excess energy over threshold in our case, i.e. 
5M = — (rriN + 'tun*), with the corresponding value in the deuteron-like 

six-quark system J^=1+,T=0 iViV-wave), = M^e — 2mN- In the latter 

case, as we know, there exists no compact 6q-state and the bound system (deuteron) 
consists of two weakly bound nucleons, which with overwhelming probability are far 
from each other. 

We see from table ^ that 6M is very similar to 5Ms&. This maybe a hint, that 
there could be a weakly bound state (below the N(939)+N*(1535) threshold) in the 
J^=0“,T=0 system. The reason is quite clear. There are more symmetrical pairs 
with CS-symmetry, [321]c5', in the J'^=0“,T=0 system, than with the [2^]cs symme¬ 
try corresponding to the deuteron-like s® conhguration. But the ST Young pattern 
for the J^=1+,T=0, s® conhguration, [ddlsr, is more symmetrical than [321]5'7-, in¬ 
herent in the state (1^) according to eq.(^). As a consequence the gain in energy 
coming from the color-magnetic forces (prefering a symmetric CS-con£guration) is 
approximately compensated by the loss in energy arising from the chiral-exchange 
interaction (which prefers symmetrical ST-con£gurations). As a consequence, both 
systems have very similar excess energies. 

In table ^ we have shown the excess 5M calculated with the same harmonic 
oscillator parameter feg in fhe six-quark system as for the baryons fog = ^n- How¬ 
ever, with hxed parameters of the quark-quark interactions of eq.(pTD, one should 
minimize the ’’mass” eq. (0) with respect to the harmonic oscillator parameter feg 
in the six-quark trial wave function. The corresponding value Bq and the mass of 
the 6q-state are shown in the 3'^'^ and 4*'^ column of table ^ A more detailed 


discussion of these results for different parameter sets can be found in 
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5 Configuration mixing results 


Here, we report on the diagonalisation of the Hamiltonian (p!6[) in the space includ¬ 
ing 11 harmonic oscillator configurations (the only state in the N=1 shell, and 10 
states with a total number of three harmonic oscillator quanta excited: N=3,[42]x, 
L,S,T=1,1,0). The quantum numbers and our labelling of these states, as used in 
the calculation, are given in table 0. All the needed matrix elements 


H,, = {i\H\j) 




(41) 


can be found in appendix B. 

It is important to mention the well known fact that a h-type attractive two-body 
potential would result in a collapse in the three-body system, if a complete basis 
were used. Thus, the h-type forces in eq.(|l5|) have to be considered as an effective 
interaction, only valid for a given finite basis. We included the contact forces in 
perturbation theory and in the full diagonalisation, and obtained within 2 to 3 MeV 
the same masses for the calculated resonance J^=0“,T=0. 

Since our basis is not complete, we improve its flexibility by using the harmonic 
oscillator parameter in the six-quark basis functions as a nonlinear variational 
parameter. More clearly, with fixed parameters of the quark-quark interactions, we 
minimize the lowest eigenvalue of (^T]) with respect to b^. This means in particular, 
that the harmonic oscillator parameters (and thus the root mean square radius) of 
the baryons b^ and in the dibaryon fcg need not to be the same. 

Let us now discuss different contributions to the mass from the Hamiltonian 


(pT|). Since the diagonalisation of eq. (^TD is a very nonlinear procedure, it is dif- 
hcult to estimate the importances of the separate parts of the Hamiltonian ( 0 )- 
Nevertheless, diagonalising for example the conhnement matrix (table |B.1|) alone, 
one recognizes that the mixing due to the confinement is essential. This result is 
independent of the chosen parameter set, since the parameter dependence factors 
out for the diagonalisation as is seen from eq.(p.3|) in appendix B. 
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The same argument holds for the individual diagonalisation of the color-Coulomb 


interaction eq. (|B.4|) . In this context, it is important to see, that large mixing ampli¬ 
tudes for the diagonalised system do not necessarily lead to a large gain in energy 
by the mixing calculation. For the conhnement for example, the energy gain is 
rather important (around 100-150 MeV) whereas in the case of the color-Coulomb 
interaction, we gain only about 30 MeV. 

One should mention here, that the effect of all forces, that involve only spatial 
and color degrees of freedom in their potentials, could be studied more easily in the 
basis with intermediate S'T-symmetry. Here, these forces are diagonal with respect 
to the spin-isospin symmetry, and we have no mixing between states with different 
S'T-symmetries for this part of the Hamiltonian. This means that only two of the 10 
chosen states from the N=3 shell with [321]5'r, (A/i) = (11) and [321]5'r, (A/r) = (30) 
would be mixed in a ST-basis to the iV = 1, [321]5r state by the confinement, 
color-Coulomb and cr-exchange forces. Using the ST reduction chain (PD|), only 
the color-magnetic, tensor and/or pion interactions would mix among states with 
different ST. 


With all parts of the Hamiltonian (hOD , this analysis of the interplay of the 
different potential parts cannot be done anymore. Of course, the difference in kinetic 
energy from the N=1 and the N=3 shell of ~300 - 500 MeV is the most obvious 
reason for the rather moderate gain in energy as seen from table ^ table shows 
in its last 2 columns hrst the harmonic oscillator length b^, for which the minimum 
of the mass is reached, and second the corresponding mass eigenvalue M for the 
chosen parameter set. This value should be compared (for the parameter sets I to 
IV) to the N(939)-|-N*(1535) threshold of 2474 MeV. 

Comparing the results with and without conhguration mixing, we see that the 
inclusion of the excited states shifts the mass down by about 100 - 150 MeV. In 
the best case (for set H), the mixing result is 90 MeV above the N(939)N*(1535) 
threshold. The inclusion of higher conhgurations makes the 6q-system larger, and 
the conhnement strength determines essentially the size of the dibaryon h^. There- 
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fore, since parameter set V allows for a large dibaryon of fee = 1-25 fm, and since the 
kinetic energy is proportional l/&^, we see qualitatively, why this set gives the lowest 
’’mass” of the dibaryon in our basis. One can see in table that this parameter set 
does not allow for a description of the excited baryon N*(1535). Thus, we conclude 
that it is difficult to describe the baryon spectrum and the dibaryon mass with the 
same confinement parameters. 

6 The J^=0 ^ T—2 SIX cjUcirk st3.tG cis ci CciiididcitG 
for thG d’ 


In this section, we discuss shortly the question of a J'^=0“, T=2 six-quark state, 
since from the very beginning, it is not clear which state, T=0 or T=2, is lower in 
mass 1^. Both possible states can not decay into NN-channels due to the Pauli- 
principle. 


In our approach, following the classification in section 3, we get two states 
with one excitation quantum N=1 and the quantum numbers J^=0“, T=2. 


|1) = |iV=l,[51]x,(10),L = l,^=l,T = 2,[321]xs) 

|2) = \N=l,[51]x,m,L = l,S = l,T = 2,[42]xs) 


(42) 

(43) 


With the other possible reduction chains (p9|) and (^), we get for (pOD the two 
states [42]5'r and [321]5'r, and for (|^) [31^]c5 and [21‘^]c’5. The necessary diagonal 
matrix-elements for the states (42) and (|43|) are given in appendix C. 

In columns 5 and 6 of table we give the minimized expectation values (|C.1|) 
and (|C.2|) for the lowest possible states |1) and |2) in comparison with the result 
for the N=l, J^=0“, T=0 state. From table ^ we see, that within the constituent 
quark model both the J^=0“,T=0 and the J'^=0“,T=2 state are nearly degenerate. 

The reason for this degeneracy is clearly the competition between the color- 
magnetic part of the gluon-exchange and the central pion-exchange contribution. 
In analogy to the discussion in section 3, we see here, that the T=0 state with 
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its [321]cs symmetry provides more attraction from the color-magnetic interaction, 
than the less symmetric C'S'-symmetries of the T=2 states. The competing process 
is the pion-exchange part, which provides for the [42] 5 ^ symmetry of the T=2 state 
more gain in energy, than in the N=l, T=0, [321]5'r state. So, at the end, the results 
depend very much on the given parameter set, or more precisely on the relative 
importance of the one-gluon- and the pion-exchange. Both six-quark systems are 
quite close in energy, and from the point of view of the constituent quark model 
neither the T=0 nor the T=2 state are energetically favoured. 


7 Conclusion 


We have investigated the six-quark system with the quantum numbers J^=0“, T=0 
within the constituent quark model with and without chiral interactions (vr- and cr- 
exchange) between quarks. The aim of this study was to look for a possible dibaryon 
with these quantum numbers. It was argued [j^, that such a dibaryon is seen as a 
narrow peak in double charge exchange reactions ( 7 r''', 7 r“) on various nuclei. 

We have classified all possible six-quark states J'^=0“, T=0 with N=1 and N=3 
harmonic oscillator excitations using different reduction chains. We then have diag¬ 
onalized the microscopic Hamiltonian including the linear quark-quark confinement, 
chiral interactions between quarks, and the effective one-gluon exchange-potential 
in the basis consisting of the lowest (N=l) state and the 10 presumably most im¬ 
portant conhgurations from the N=3 shell. We have argued, that the calculated 
mass-eigenvalue should be compared to the 2-baryon N(939)-|-N*(1535) threshold. 

We have found, that if one fits all parameters to describe the nucleon and its 
lowest resonances, the mass of such a ’’dibaryon” lies above the N(939)-|-N*(1535) 
threshold, when taking our restricted basis of 11 states. Since we gain around 100 
~ 150 MeV by admixing the 10 chosen states, it seems not to be unrealistic, that 
a more complete basis would lead to a ’’dibaryon mass” below the N(939)N*(1535) 
threshold of 2474 MeV. So, our calculations do not deny the existence of a dibaryon 
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with J^=0“, T=0. Its mass could be some 10 MeV’s below the NN* threshold 
(within the framework of our model). But the discrepancy with the ’’experimentally 
found” mass is still large. 

We cannot claim for sure, that the constituent quark model ’’excludes” the in¬ 
terpretation of the peak in (7r+,7r“) reactions as a signal from the dibaryon with 
the above mentionned quantum numbers, since our basis is not complete, and as a 
consequence, we lose probably some of the possible short-range and few-quark cor¬ 
relations in the six-quark system. Another point to mention is the following: The 
6q-system differs qualitatively from the description of a baryon in the sense, that 
the 6q-color symmetry is a mixed symmetry, in contrast to the fully antisymmet¬ 
ric color-wavefunction of a baryon. Probably, the conhnement for a 6q-system (in 
addition to nuclear medium effects in the nucleus) differs qualitatively from the con¬ 
finement mechanism in the baryons. At this point, we should add, that the choice 
of a quadratic or a so-called color-screened exponential 


rather than a linear 


conhnement potential does not change qualitatively our results. If the assumption 
of a dibaryon resonance at 2065 MeV is conhrmed, we could probably improve our 
understanding of conhnement in multi-quark systems. 

We have also analyzed a possible J^=0“, T=2 dibaryon. The corresponding 
lowest shell model state lies below the J^=0“, T=0 one, provided that the one- 
pion-exchange forces play an essential role in the 6q-system. 

Finally, we have shown, that we could describe an ’’observed dibaryon” with a 
mass close to Md'= 2065 MeV, by reducing the conhnement strength Oc (string ten¬ 
sion), down to Oc — 25 MeV/fm [|^ (see parameter set V). But, such a conhnement 
strength does not give the correct energy for the negative parity one-quantum exci¬ 
tation N*(1535). Thus, we observe two opposing tendencies within our model: We 
need a small conhnement strength (and therefore a large size parameter for the 
nucleon) to get a rather large dibaryon, which is consequently rather light in mass. 
On the other hand, we need a rather small nucleon size to describe the baryon 
spectrum, especially the N* excited states. This result suggests, that this six-quark 
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system is strongly influenced by the confinement mechanism. The above discrepancy 
can not be eliminated by another choice of radial dependence of the conhnement 
potential (quadratic or ’’color-screened” exponential) alone. If the d’ dibaryon is 
conhrmed, we suggest that the color-structure A“ ■ A“ of the conhnement potential 
should be modihed in multi-quark systems like the d’. We should stress here once 
more, that it seems impossible to describe the nucleon resonances (the 3q-sector) 
and this dibaryon resonance within the same model with the same parameters. 
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Appendix A 

In this appendix, we give the mass formulae for the baryon systems 

N(939), A(1232), and N*(1535), N**(1440) (A.l) 


within the constituent quark model, as described in the parameter htting section. 
The following harmonic oscillator basis states 


|Ao> = \N = 0,{Xfi) = {00),[3]x,L = 0,S=^,T = ^[3]sT) , 

|Ao> = |A = 0,(A/r) = (00),[3]x,T = 0,^=|T = ^[3]5T), 

Nt^> = \N = 1 ,{\ix) = {10),[21]x,L = 1,S=^,T = ^[21]st) , 
Nt^> = \N = 1,{\^) = {10),[21 ]x,L = 1,S=^,T = ^[21]st) , 

|Ai> = I A = 2,(A/i) = (20),[3]x,T = 0,^=^,T = ^[3]5T), 

\N2> = I A = 2,(A/i) = (20),[21]x,T = 0,^=i,T = l[21]5T), (A.2) 
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are needed for the calculation of the appropriate matrix elements which determine 
the masses of the baryons N(939), A(1232), N*(1535) and N**(1440) in our chosen 
parameter htting scheme. The ground states are denoted by the subscript ” 0”. The 
detailed wave functions can be found in the second paper of [^]. 

Our Ansatzes for the different baryons are: 


\N) = |7Vo) : |A) = |Ao) 

|jv) = A ly-')+i^vh) 

\N") = a\Nj)+l3\N2} (A.3) 


To be rather short, we give here all necessary (diagonal and off-diagonal) matrix 


elements in their analytical form in one table, 
h,/ = Snigdij + t 




mqb'^ 

. , a 

+ [ceij + cmij) 


-|- dij —— — 2C/ Oij “h d 


CXc 


y/^rrigb^ 


+ 14J + + VS 


+ 


(A.4) 


The notations for the pion central-contribution are the same as in ([B.9|) and for the 
cr-exchange the notations are the same as in ( p.l0|) . As in appendix B, we do not 


give the lengthy formulae for the pion-tensor contribution in explicit form. 



Appendix B 

In this appendix, we present all matrix elements 

Hji = Hij = {i\H\j) ; i, j = 1... 11 

= Qrrig 5ij + + 1/" + 1/").. (B.l) 
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for the different parts of the Hamiltonian eg. (p!6D . in the tables |B.1| throngh |B.4| , 
excluding the pion tensor contribution again. As mentioned above, the labelling of 
the basis states, as they are used here, is given in table 


The kinetic energy contribution is given by 

^ ii i = j = 1 


rpkin 


4 m„ 6^ 
21 1 
4 mq 62 


if i = j = 2 ... 11 
if i^j 


(B.2) 


The coefficient Uij of the conhnement contribution 

vr' = {*1 E (a.r - C) \j} = a., ^ - 4C4 

k<l ^ 

is given by 






(B.3) 


Tables |B.2a| - [B.2d| give the different parts of the one-gluon exchange potential, 
starting with the color Coulomb contribution 


yCC _ /n V — • A“ hi) - h 


k<l 


The next table gives the so-called color-electric part 


k<l ^ '"‘q 


da 


ce. 


V 


This is the color-magnetic contribution 

^ X—> cr.s 71 

The next table gives the gluon tensor contributions 




Oia 


6 


k<l 


Xl • A? Ski h-\j) = cU, - 

Ibm^ 


Introducing the following notations 

poo rp‘^ 

I(n){jn^,bN) = r’"exp(-^ - m^r)dr 


/(o)(m^,6jv) = ftjv erfc(^2^^^ exp( 

2 

erfc( 2 ;) = —= / exp{—x^)dx 

a/tT J z 


rnSSP' 




d 


gn+l) 


dm^. 


'(n) 
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for the radial integrals needed in the calculation of the Yukawa-like pion- and 
sigma-potential matrix elements, allows us to present the chiral contributions in 
the following form: First, for the pion central contribution, we write 


A2 


j^2 _ ^2 4^ V vr 


9 P^J ./^ r, , , .Q Y3) , . H5) \ ,2 




(B.9) 


So, in the four next tables (table |B.3a) - [B.3d| ), we give the coefficients p,pl,p3,p5 
respectively. Let us mention once more, that since the pion tensor contributions in¬ 
volve all above integrals from J(i) up to /(s), and gives rather negligible contributions, 
we do not give the analytical expressions here. 


Finally, table p.4| concerns the a-potential, where we give only the non-zero 
matrix elements. The notations are the same as for the pion, i.e. 

A^ Sn l2 (/ _ ^ „ 1(1 
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Appendix C 


In this appendix, we present the expectation values of the Hamiltonian (hGl) for 
the two states of eq.(0^ and eq.( 
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One recognizes the CC-term ~ ^ and the added CE and CM-term ~ The 
second ^-term comes from the tensor part of the one-gluon exchange potential. 
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With the notations of appendix B, eq. ([B.9|) and eq.( |B.10|) , the central parts of the 
chiral interactions are given in the following list 


VSH) 

1/^(1) 

^c(2) 

V^{2) 
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The tensor contributions of the pion are given by 
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Pt(i) = ^ 


Pt(2) = 
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i>N 

A 

nig 


dc 

C 

iV*(1535) = XN[~^ + 

]\f** 

Set 

[fm] 

[fm“^] 

[MeV] 


[MeV/fm] 

[MeV] 

A 

A 

V*(1535) 

V*(1650) 

[MeV] 

I 

.45 

5.07 

338 

.127 

462 

549 

.98 

.22 

1535 

1697 

1678 

II 

.47 

7.6 

296 

.074 

423 

500 

-.96 

.30 

1535 

1713 

1660 

III 

.5 

3.55 

274 

.334 

654 

676 

.998 

.06 

1535 

1683 

1797 

IV 

.5 

/ 

230 

.474 

1223 

1247 

-.86 

.51 

1535 

1844 

2120 

V 

.6 

10.14 

313 

.816 

25 

-32 

.98 

.22 

1275 

1436 

1292 


Table 1: Different sets of parameters for the quark-quark interactions of eqs.(10)-(15), fitted as 
described in the text. The hrst 6 columns show the harmonic oscillator length used to describe 
the baryons, the cutoff A parametrizing the hnite size of the pion-quark and sigma-quark vertex, 
the constituent quark mass rUg, the strong coupling constant o;* of the one-gluon exchange and 
the strength (slope) Oc and offset C of the conhnement potential, respectively. The diagonalisation 
of the mass matrix in the basis of two states n[ ^ and N 2 ^ gives the energies of the two negative 
parity resonances N*(1535) and N*(1650). The next four columns give the amplitudes and masses 
of these two resonances, respectively. The last column gives for illustration the mass of the lowest 
two-quantum excitation (N=2) for the parameters given, calculated in a basis of the two most 
important states. 


34 






















































































[f]x 

|51| 

|51| 

[42] 

[42] 

[42] 

[411] 

[411] 

[321] 

|321] 

[3111] 

(A/i) 

(11), (30) 

(11) 

(11), (30) 

(11) 

(30) 

(11), (30) 

(11) 

(11) 

(11) 

(00) 

LST 

110 

220 

110 

220 

330 

110 

220 

110 

220 

000 

[f]cs 

|321] 

13211, 

pn^i 

[42], [321], 
[31^], [23], 
[211 

|321] 

[23] 

[321], 

PC] 

[321], 

|2“e] 

[42],[321]i, 

[321]2,[313], 

[211 

[321]i, 

[321]2, 

[2212] 

[412] 

[f]sT 

|321] 

[42], 

|321] 

[51], [411], 
[33], [321], 
[2212] 

|321] 

133] 

[411], 

1321] 

[42], 

[321] 

[51], [411], 
[321]i,[321]2, 
[2212] 

[42], 

[321]i, 

[321]2 

[311 

[f]xs 

[42] 

[6], 

|42] 

[6],[42]i, 

[42]2,[313], 

[222] 

[42] 

PI 

[42], 

[313] 

[42], 

[31=] 

[42]i,[42]2, 

[3l1i,[3ll2, 

[222] 

[42], 

[311, 

[222] 

[311 

states 

2 

2 

10 

1 

1 

4 

2 

5 

3 

1 


Table 2: Here, we give a complete list of all possible six-quark states with three internal harmonic 
oscillator quanta N=3 and the quantum numbers J'^=0“, T=0. The hrst line dehnes the spatial 
permutational symmetry [f]x, and the second the corresponding possible Elliot symbols (A/i) for the 
harmonic oscillator. The third line shows the total orbital angular momentum L, the total spin S and 
the total isospin T. The next three lines give the intermediate permutational symmetry needed for an 
unambigious classihcation of the states within a certain reduction scheme. Line 4 corresponds to the 
reduction chain (29) with explicit intermediate color-spin symmetry, line 5 to (30) with intermediate 
spin-isospin symmetry and line 6 to (31) with intermediate spin-orbital symmetry. The redaction 
chain (31) was nsed in the present calculation. The last line gives the nnmber of states in each colnmn, 
respectively. 


35 












































































excess energies 

groundstate 

T=2 

mixing results 

■ 

6M 

{Bq = Bn) 

{Be = Bn) 

be 

/\/f'r=o 

^yJ-N=i 

II 

to 


be 

M 

Set 

[MeV] 

[MeV] 

[fm] 

[MeV] 

[MeV] 

[MeV] 

[fm] 

[MeV] 

I 

525 

509 

.59 

2705 

2564 

2744 

.60 

2588 

II 

566 

555 

.65 

2680 

2530 

2705 

.65 

2565 

III 

624 

603 

.60 

2898 

2812 

2949 

.62 

2760 

IV 

850 

924 

.56 

3288 

3303 

3333 

.59 

3114 

V 

441 

400 

1.24 

2162 

2146 

2169 

1.24 

2132 


Table 3: For different sets of parameters I to V defined in table 1, the first column 
5M gives the mass difference of the calculated lowest J^=0“, T=0 state and the 
N(939)N*(1535) threshold. For comparison, the second column 5Ms& gives the mass 
excess for six quarks, coupled to deuteron-like quantum numbers, in one harmonic 
oscillator s® with respect to the N(939)N(939) threshold at 1878 MeV. Both mass 
excesses in columns 1 and 2 are calculated for the oscillator length equal to 
the value hM minimizing the nucleon mass. Columns 3 and 4 show the oscillator 
length obtained by minimizing the single N=l, J^=0“, T=0, [51]x, (A/r) = (10), 
[f']xs = [42] conhguration and the corresponding mass. Columns 5 and 6 list the 
masses of the two N=l, J'^=0“, T=2 states. Here again, the oscillator length (not 
given) is varied to minimize both energies separately. Columns 7 and 8 give the 
hnal result for the oscillator length Bq, minimizing the lowest J^=0“, T=0 energy 
eigenvalue, and the corresponding mass M. This mass should be compared with the 
observed experimental resonance energy of 2065 MeV. 


harm.osc.j) 

i = 1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

N 

1 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

[f]x 

|51] 

[42] 

[42] 

[42] 

[42] 

[42] 

[42] 

[42] 

[42] 

[42] 

[42] 

(A/i) 

(10) 

(30) 

(11) 

(30) 

(11) 

(30) 

(11) 

(30) 

(11) 

(30) 

(11) 

L, S, T 

1, 1, 0 

[f']xs 

[42] [6] [6] [2®] [2®] [31®] [31®] [42]^ [42]i [42]^ [42]^ 


Table 4: Quantum numbers and our labelling of the states included in the conhgura¬ 
tion mixing to calculate the lowest J'^=0“, T=0 six-quark conhguration. The hrst 
line gives our labelling of the states. The next line lists the respective number of 
excited harmonic oscillator quanta N, and line three represents the orbital permuta- 
tional symmetry of the Young tableaux. The fourth line shows the Elliot quantum 
numbers of the harmonic oscillator (A/r), while line hve gives the total orbital an¬ 
gular momentum L, the total spin S and the total isospin T of the conhgurations. 
The last line denotes the intermediate symmetries of the Young tableaux [f']xs in 
the orbital-spin space. 
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Table B.l: Analytic expressions for the confinement matrix elements 



Table B.2a: Analytic expressions for the color coulomb matrix elements 
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Table B.2c: Analytic expressions for the color magnetic matrix elements ((5-force) 
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Table B.2d: Analytic expressions for the gluon tensor matrix elements 
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Table B.3a: Analytic expressions for the pion exchange matrix elements (central part) 
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Table B.3c: p3 for the analytic expressions for the pion exchange matrix elements 
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Table B.3d: p5 for the analytic expressions for the pion exchange matrix elements 


(ij) 


Slij 

s3ij 

s3ij 

(11) 

1 

12 

1 

0 

(22), (44), (66), (88), (10,10) 

IH 

29 

H 

mim 

(33), (55), (77), (99), (11,11) 

IH 

28 

H 

H 

(23), (45), (67), (89), (10,11) 

lEE 

1 

3 

I 

15 


Table B.4: Analytic expressions for the non-zero sigma exchange matrix elements 
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